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Abstract 

The XOR-Ising model on a graph consists of random spin configurations on vertices of 
the graph obtained by taking the product at each vertex of the spins of two independent 
Ising models. In this paper, we explicitly relate loop configurations of the XOR-Ising 
model and those of a dimer model living on a decorated, bipartite version of the Ising 
graph. This result is proved for graphs embedded in compact surfaces of genus g. Using 
this fact, we then prove that XOR-Ising loops have the same law as level lines of the 
height function of this bipartite dimer model. At criticality, the height function is known 
to converge weakly in distribution to a Gaussian free field |dT07bj . As a consequence, 
results of this paper shed a light on the occurrence of the Gaussian free field in the 
XOR-Ising model. In particular, they provide a step forward in the solution of Wilson's 
conjecture [Willi] , stating that the scaling limit of XOR-Ising loops are level lines of the 
Gaussian free field. 



1 Introduction 

The double Ising model consists of two Ising models, living on the same graph. It is related 
[KW7H fWu7H [Fan721 Weg72| to other models of statistical mechanics, as the 8- vertex model 



[Sut70l IFW70] and the Ashkin- Teller model [AT43| . In general, the two models maybe 
interacting. However, in this paper, we consider the case of two non-interacting Ising models, 
denned on the dual G* = (V* ,E*) of a graph G = (V,E), having the same coupling constants 
(J e *) e * g £*, where the graph G is embedded either in a compact, orientable surface E of genus 
g > 0, or in the plane. 

We are interested in the polarization of the model [KB79], also referred to as the XOR- 
Ising model [Willi] by Wilson. It is denned as follows: given a pair of spin configurations 
(a, a') G {—1,1}^ x {— 1,1}^ , the XOR-spm configuration belongs to {—1,1}^ and is 
obtained by taking, at every vertex, the product of the spins. The interface between ±1 
spin configurations of the XOR-configuration is a loop configuration of the graph G. Using 
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extensive simulations, Wilson [Willi] finds that, when G is a specific simply connected domain 
of the plane, and when both Ising models are critical, XOR loop configurations are contour 
lines of the Gaussian free field, with heights of the contours spaced y/2 times as far apart 
as they should be for the double dimer model on the square lattice. Similar conjectures, 
involving SLE rather than the Gaussian free field, are obtained through conformal field theory 
[IR11[ iPSllj . Results of this paper give a partial proof of the occurence of the Gaussian free 
field. 

The first part of the paper concentrates on finite graphs embedded in surfaces. We explicitly 
relate XOR loop configurations to loop configurations in a bipartite dimer model, implying 
in particular that both loop configurations have the same probability distribution. In the 
second part, we prove that this correspondence still holds for a large class of infinite planar 
graphs, the so-called isoradial graphs [Ken02, KSQ5], at criticality, and make the connection 
with Wilson's conjecture. Here is an outline. 

Outline 

Section [2j One of the tools required is a version of Kramers and Wannier's low/high 
temperature duality K\V 1 lal IKW41b| in the case of graphs embedded in surfaces of genus 
g, with boundary. In the literature, we did find versions of this duality for graphs embedded 
in surfaces of genus g [L G94] . but we could not find versions taking into account boundaries. 
This is the subject of Propositions [6] and [7J it involves relative homology theory and the 
Poincare-Lefschetz duality. 

Sections [3] and [4] consist in the extension to general graph embedded in a surface of genus 
g of an expansion due to Nienhuis [Nie.84] , which can be summarized as follows. Consider 
the low temperature expansion of the double Ising model, i.e. consider pairs of polygon 
configurations separating clusters of ±1 spins of each spin configuration. Drawing both 
polygon configurations on G yields an edge configuration consisting of monochromatic edges, 
that is edges covered by exactly one of the two polygon configurations, and bichromatic edges, 
that is edges covered by both polygon configurations. Monochromatic edge configurations 
exactly correspond to XOR loop configurations, and separate the surface S into connected 
components £i,--- , Sjv- Inside each connected component, the law of bichromatic edge 
configurations is that of the low temperature expansion of an Ising model with coupling 
constants that are doubled. As a consequence, the partition function of the double Ising 
model can be rewritten using XOR loop configurations and bichromatic edge configurations, 
see Proposition [TUl 

Fixing a monochromatic edge configuration, and applying low/high temperature duality to 
the single Ising model corresponding to bichromatic edges, yields a rewriting of the double 
Ising partition function, as a sum over pairs of non intersecting polygon configurations of 
the primal and dual graph, where primal polygon configurations exactly correspond to XOR 
loop configurations, see Proposition [T2] and Corollary [T31 Note that there are quite a few 
difficulties in the proofs, due to the fact that we work on a surface of genus g. 
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Proposition 1. 

The double Ising partition function for a graph embedded on a surface of genus g can be 
rewritten as: 

( 2e~ 2J ' 

Zd-ism g (G* , j) =ci ^2 ( n 1 

{(P,P*)eP (G)xP (G*): VeSP 

PnP*=0} 

where primal polygon configurations of V°(G) are the XOR loop configurations, and C\ = 
2 |y * l+2s+1 (n e6 pCo S h(2J e ,)). 

Section O In Section we define the 6- vertex model on the medial graph G M constructed 
from G. Reformulating an argument of Nienhuis [Nic84] , we prove that the 6- vertex partition 
function can be written as a sum over non-intersecting pairs of polygon configurations of the 
primal and dual graph. 

In Section 15. 2\ we define the dimer model on the decorated, bipartite graph G^ constructed 
from G. Then, we present the mapping between dimer configurations of G^ and free fermionic 
6- vertex configurations of G M [WL75} IDubllb] . Using both mappings, one assigns to every 
dimer configuration M a pair Poly(M) = (Poly 1 (M), Poly 2 (M)) of non intersecting primal 
and dual polygon configurations. The weights of the 6- vertex model chosen to match those 
of edges in the mixed contour expansion of the double Ising model satisfy the free fermionic 
condition. As a consequence, we then obtain, see also Proposition 1 17t 

Proposition 2. The dimer model partition function Z^ mer {G^, J) can be rewritten as: 

z 3lm „ (0 v) = 2 e (n if^r) ( n ItfS^ 

{(P,P*)EV°{G)xV°(G*): VeSP / \e*eP~ 

PnP*=0} 



where primal and dual polygon configurations ofV°(G) x V "(G*) are the Poly configurations. 
Combining Proposition [1] and Proposition [2] yields the following, see also Theorem 1181 



Theorem 3. XOR loop configurations of the double Ising model on G* have the same law as 
Poly 1 configurations of the corresponding dimer model on the bipartite graph G® : 

VP e V°(G), P d . Ising [XOR = P] = F° Q [Poly 1 = P] . 

Remark 4. In the paper [Dubllb] , Dubedat relates a version of the double Ising model 
and the same bipartite dimer model in two ways. The first approach uses explicit mappings, 
most of which are present in the physics literature, and goes as follows. Consider a slightly 
different version of the double Ising model, with one model living on the primal graph G 
and the second on the dual graph G* . This double Ising model can be mapped to an 8- 
vertex model [KW7H IWu71] on the medial graph. Using Fan and Wu's abelian duality, 
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this 8-vertex model |FW70] can be mapped to a second 8-vertex on the same graph. When 
coupling constants of the two Ising models satisfy Kramers and Wannier's duality, the second 
8-vertex model is in fact a free-fermionic 6-vertex model. The free-fermionic 6-vertex model 
can in turn be mapped to a bipartite dimer model, a result due to [WL75J in the case of the 
square lattice, and extended by [Dub lib) in the general lattice case. It can also be seen as 
a specific case of the mapping of the free-fermionic 8-vertex model to a non-bipartite dimer 
model of [FW70]. Note that this bipartite dimer model is the model of quadri-tilings studied 
by the second author in [dT07a] and [dT07b] . 

When performing the different steps of the mapping, Dubedat keeps track of order/disorder 
variables, in the vein of |KC71j . Using results of a previous paper of his [Dub 11 a] , this allows 
him to compute critical correlators in the plane. For simply connected regions, this result 
has independently been obtained by Chelkak, Hongler and Izyurov [CHI12J. 

Our goal here is different, since we aim at keeping track of XOR-configurations. This infor- 
mation is lost in the above approach. Indeed Fan and Wu's abelian duality for the 8-vertex 
model can be compared to a high temperature expansion, where configurations cannot be 
interpreted using the initial model. 

The second approach uses transformation on matrices. The partition function of the double 
Ising model can be expressed using the determinant of the Kasteleyn matrix of the Fisher 
graph |Fis66j ; whereas the partition function of the bipartite dimer model can be expressed 
using the Kasteleyn matrix of the graph . Dubedat shows that the two matrices are related 
through transformations not affecting the determinant. Using the fact that the partition 
function of the double Ising model is also related to the determinant of the Kac-Ward matrix 
|KW52j . Cimasoni and Duminil-Copin use the same approach to relate the Kac-Ward matrix 
to the matrix of the same bipartite dimer model [CD12]. Their purpose is to identify the 
critical point of general bi-periodic Ising models, see also Li [LilO] ILil2j for the case of the 
square lattice with arbitrary fundamental domain. 

When performing the above transformation on matrices, one again looses track of double 
Ising configurations, and in particular of XOR-configurations. 

Using Nienhuis' mapping [Nie84|, the main contribution of this paper is to provide a coupling 
between the double Ising model and the bipartite dimer model, which keeps track of XOR 
loop configurations, and is valid for graphs embedded in surfaces of genus g. 

Section [6j Suppose now that the two Ising models are critical and defined on the dual of 
an infinite isoradial graph G, see Section 16.11 for definitions. Then, the dimer model on the 
corresponding graph G® is also critical in the dimer sense. Using the locality property of 
both probability measures on Ising [ BDTllj . and dimer configurations [dT07bj on isoradial 
graphs at criticality, we prove that the equality in law stated in Theorem [3] still holds in this 
infinite context. See Theorem 1201 

Section The graph G® being bipartite, using a height function denoted h, dimer configu- 
rations can naturally be interpreted as discrete random interfaces. Our second theorem, see 
also Theorem 123} proves the following 
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Theorem 5. XOR loop configurations of the double Ising model defined on G* have the same 
law as level lines of the restriction of the height function h to vertices of the dual graph G* . 

Theorem [5] can be seen as a proof of Wilson's conjecture mentioned above (see 17.21 for a 
precise statement) in the discrete setting, since it is known that the height function h, seen 
as a random distribution, converges in law in the scaling limit to times the Gaussian free 
field in the plane |dT07bj . In particular, we explain in this section the special value of the 
spacing. It would yield a complete proof of the conjecture if we could overcome the same 
technical obstacles for the proof of the convergence of double dimer loops to CLE4. 

Acknowledgments: We would like to thank warmly Thierry Levy for very helpful discussions 
on relative homology. 



2 Ising model on graphs embedded in surfaces with boundary 

In the whole of this section, we let G be a graph embedded in a compact orientable surface 
of genus g (g > 0), and G* be its dual graph. The embedding of G* is chosen such that dual 
vertices are in the interior of the corresponding faces. 

Fix some positive integer p, and for every i £ {0, • • • ,p — 1}, let Bi be a union of closed faces 
of G homeomorphic to a disc, such that nf=o &i = Denote by E the surface of genus g 
from which the union of the IV s is removed. Then E is a compact orientable surface of genus 
g, with boundary dT, = OBq U • • • dB p -\. 

Let Gy = (Vs, Ey) be the subgraph of G defined as follows: Vy consists of vertices of V n E; 
and Ey consists of edges of E joining vertices of Vy, from which edges on the boundary dT 
are removed. Let G^ = (V£,E^) be the subgraph of G* whose vertices are vertices of V* n E, 
and whose edges are edges of G* joining vertices of V£, see Figure [1] for an example. Note 
that the graph G^ contains all edges dual to edges of Gy, i-e. there is a bijection between 
primal edges of Gy and dual edges of G^ . 

Fix a collection of positive constants ( J e * ) e *eE* attached to edges of G* , referred to as coupling 
constants. The Ising model on G s with coupling constants (J e *) is defined as follows. A spin 
configuration a of is a function of the vertices of V^j with values in {— 1,+1}. The 
probability of occurrence of a spin configuration a is given by the Ising Boltzmann measure, 
denoted Pisingj and defined by: 



G {-1, Pi sing (a) = / exp [ £ J e *a u *a v * 

Rising V^S' J ) 



where Zisi ng (G s , J) = exp I J e *o~ u *o~ v * J is the Ising partition function. 

<76{-l,+l} V S \e*=u*v*&E^ J 

Note that to simplify notations, the inverse temperature is included in the coupling constants. 
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Figure 1: The graph G is a piece of 1? embedded in the torus. The union of faces (-Bi)«e{o.i,2} 
is pictured in light grey. The graph Gs consists of black plain lines, and the dual graph 
of black dotted lines. The paths (7i)f =1 defining defects of Section 12.11 are drawn in thick 
black lines. 



2.1 Ising model with defect lines 

Let us denote ./V = 2g + p — 1. We now introduce 2^ versions of the original Ising model. Let 
7ii • • • ;7iV be N unoriented paths consisting of edges of the primal graph (?£, see Figure [1] 
for an example, where 

• for i G {1, ... ,g} , the paths 72i-i,72« wind around the i-th handle in two transverse 
directions, 

• for i G {1, . . . ,p — 1}, the path ~f2g+i joins OBq and dBi. 

Let e be one of the 2^ possible unions of paths {J^jfi, where I C {1, • • • , N}. Then, we 
change the sign of coupling constants of dual edges crossing paths of e. Spin configurations 
are defined as above, and so is the probability measure on spin configurations. This defines 
the Ising model with coupling constants (J e *) and defect condition e. 

In fact, the appropriate framework for defining the Ising model with defects, is relative 
homology theory, see Appendices IA.2IA.3I and IA.4I The first homology group o/S relative to 
its boundary <9E is denoted by dT,; Z/2Z). The collection of paths (71, • • • , 7at) defined 

above, is a representative of a basis T = (71,- •• , 7jv) of the first relative homology group 
fli(£,0£;Z/2Z). 

Let e denote the relative homology class of e in dT,; Z/2Z). Then, it will be clear from 

the low temperature expansion of Section 12.21 that the partition function is independent of 
the choice of basis and of the choice of representative of e. As a consequence, we refer to this 



6 



model as the Ising model with coupling constants (J e *) and defect condition e, and denote by 
Zf sing (G^, J) the corresponding partition function. Nevertheless, since we want the change 
of signs of coupling constants to be well defined throughout the paper, we fix representatives 
of relative homology classes in Hi(E, dT; Z/2Z), using the collection of paths 7i,--- ,7at 
defined above. Note that the original Ising model introduced has defect condition e = and 
e is the empty collection of paths. 

2.2 Low and high temperature expansion 

Proposition[U]below extends the low temperature expansion of Kramers and Wannier |KW41a[ 
IKW41b| to the case of graphs embedded on a compact orientable surface with boundary. It 
consists in rewriting the Ising partition function as a sum over polygon configurations of the 
graph Ge, 'separating' clusters of ±1 spins, see Figure [2] (left) for an example. 




Figure 2: Left: polygon configuration of Gy corresponding to a spin configuration of the 
Ising model with defect condition e = 74. Right: polygon configuration of G^. 

A polygon configuration of Gy is a subset of edges of Gy, such that vertices not on the 
boundary dT are incident to an even number of edges. There is no restriction for vertices on 
the boundary <9X. Let us denote by V(Gy) the set of polygon configurations of Gy- 

Let e be a generic element of Hi(E, dT; Z/2Z), and let V e (Gy) denote the set of polygon 
configurations of Gy whose relative homology class in Hi(T, dT; Z/2Z) is e. This defines a 
partition of 'P(Gs): 

V(Gy) = |J V e {GY). 

ee-Hi(£,dE;Z/2Z) 

Proposition 6 (Low temperature expansion). 

For every relative homology class e E e?E;Z/2Z), 

w G £< j )= 2 (n e j ^) £ n e_2Je *- w 
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Proof. Suppose for the moment that e is the class G Hi(E, <9£; Z/2Z), so that we deal with 
the usual Ising model. Using the identity ([2]) below, one can rewrite the partition function 
as a statistical sum over polygon configurations separating clusters of ±1 spins: if a u * and 
a v * are two neighboring spins of an edge e* = u*v* , then 

e J e * w = e J e * ((W= ^ } + e-^5 WutM ) . (2) 

When injecting the right hand side in the expression of the Ising partition function, the 
product over dual edges e* of e Je * can be factored out. Since primal and dual edges are 
in bijection, this can also be written as a product over primal edges. Then, expanding the 
product, we get a product of contributions for all edges separating two neighboring spins 
with opposite signs. These edges form a polygon configuration P° of separating clusters 
of ±1 spins. As a consequence P° has homology class 0, i.e. P° belongs to 

Conversely, any polygon configuration of P (Gs) is the boundary of exactly two spin config- 
urations, one obtained from the other by negating all spins, which explains the factor 2 on 
the right hand side of (P). 

Suppose now that e / 0. In the Ising model with defect condition e, coupling constants of 
edges crossing paths of the representative e are negated. For these edges, the relation ([2]) 
should be replaced by the following: 

-J e *<T u *(T v * J e * /r , -2J e * c \ 

Note that, when comparing to ([2]), the two Kronecker symbols have been exchanged. As a 
consequence, the construction of polygon configurations as above is slightly modified: the 
edge configuration, denoted by P, constructed from a spin configuration is obtained from P° 
by switching the state of every edge e in the collection of paths e, see Figure [2] (left). Then, 
the relative homology class of P in H\(Yj, <9X; Z/2Z) is: 

[P] = [J*] + [ € ] = + e = e. 

As a consequence P belongs to V e {Gy) and this, independently of the choice of representative 
of e. Conversely, any element of T >t (G's) is obtained twice in this way. □ 

For the sequel, it is useful to introduce a symbol for the sum over polygon configurations of 
the low temperature expansion. For e G Hi(E, 5S; Z/2Z), define 

zt T (Gx,J)= £ (IK 2J 4 

PeP E (G E ) eeP 

The partition function of the Ising model with defect condition e can thus be rewritten as: 

Zf sins (GhJ) = 2( J] ^*)Z£ T (G< S , J). 
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Proposition [7] below extends the high temperature expansion [K W41aj IKW41b| IWan45j to 
the case of graphs embedded in a compact orientable surface with boundary. It consists in 
rewriting the Ising partition function as a sum over polygon configurations of the graph G%, 
this time. In this case, polygon configurations do not have a simple interpretation in terms 
of spin configurations. 

A polygon configuration of is a subset of edges such that each vertex of is incident 
to an even number of edges, see Figure [5] (right) for an example. It is thus a union of closed 
cycles on G^. Let us denote by V(G^) the set of polygon configurations of G%. 

Let .Hi(£;Z/2Z) be the first homology group of S, see Appendices IA.ll IA.3I and lA.4l Then, 
to each polygon configuration of is assigned its homology class in Z/2Z). For every 

r G Hi(E; Z/2Z), we let V T {G^) denote the set of polygon configurations restricted to having 
homology class r in Z/2Z). This defines a partition of V(G^): 

V{G%) = lj V T {G^). 

T6i2i(E;Z/2Z) 

Proposition 7 (High temperature expansion). 

For every relative homology class e E Hi(E, dT,; Z/2Z), 

ZLin^Gh J) = 2™ ( n COSh(J e *)) • £ [(-1)™ Y, ( II tanh ( J - 

eeE s re#i(S;Z/2Z) P*eV T (G%.) e*eP* 

(3) 

where {r\e) is the intersection form evaluated at r and e: it is the parity of the number of 
intersections of any representative of r and any representative of e. 

For details on the intersection form, see Appendix IA.51 

Proof. This result is based on yet another way of rewriting the quantity e ±Je* (T u*< J v* f OT a 
dual edge e* = u*v* of E%. 

e ±J e *a u *a v * _ cos j 1 _j_ cr^^, s i n h J e , 

= cosh J e » (1 ± <r u * (J„* tanh J e * ) . (4) 

The partition function is expanded into a sum of monomials in (er u *) u * e y*. In the expansion, 
the spin variables come by pairs of neighbors o u *a v * and thus can be formally identified with 
the dual edge connecting u* and v* , associated with a weight ± tanh J e . . Each monomial 
is then interpreted as a subgraph of G^, the degree of a u * being the degree of u* in the 
corresponding edge configuration. Because of the symmetry a <-> —a, when resumming 
over spin configurations a, only terms having even degree in each variable remain, giving a 
factor 2 per dual vertex, and other contributions cancel. As a consequence, the contributing 
monomials correspond to even subgraphs i.e. polygon configurations of V(G^). 

We now determine the sign of dual polygon configurations. Fix r G ifi(£;Z/2Z) and a 
dual polygon configuration P* £ V T (G^). Then, edges of P* carrying a negative weight 
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are exactly those crossing paths of e. As a consequence, the sign of the contribution of P* 
corresponds to (—1) to the parity of the number of edges of P* crossing paths of e, this is 
exactly given by (r|e). The dual polygon configuration P* has thus sign (— l)( r l e ). □ 

As in the case of the low temperature expansion, it is useful to introduce a notation for the 
sum over polygon configurations of the high temperature expansion. For r € Hi(E;Z/2Z), 
define: 

ZMGhJ)= E ( II tanh(J e *)). 

p*eP T (G|,) e*eP* 

The relation between ([1]) and ([3]) can then be rewritten in the following compact form. For 
every relative homology class e G FFi(£, cffi; Z/2Z): 

Z 1t (Gx, J) = 2^*l- 1 [ [J C -^A E [(-^ZMGh J)} ■ (5) 

\ee£ s / r6Hi(E;Z/2Z) 

Using the orthogonality relation 

£ (-1)^(-1)^ = 2% T ,, 

eG-ffi(S,aS;Z/2Z) 

with 7V = 2g +p — 1, one can easily invert this relation: 

Zfi T(0i , J) . 2--l^-! + . f d£o) E [(-D W ^L T (G & J) 

\ee-B E V 7 £ePi(E,9S;Z/2Z) 

3 Double Ising model on a surface of genus g 

In the whole of this section, we let G be a graph embedded in a compact, orientable surface 
£ of genus g, and G* denote its dual graph. Since £ has no boundary, the first homology 
group Hi(S, <9£; Z/2Z) of S relative to its boundary, is identified with the first homology 
group fli(£,Z/2Z). 

Instead of one Ising model on G* , we now consider two copies of the Ising model, say a red 
one and a blue one, with the same coupling constants (J e *). These two models are not taken 
to be completely independent: we require that they have the same defect conditions, i.e. we 
ask that polygon configurations coming from the low temperature expansion of both spin 
configurations have the same homology class. 

More precisely, from the point of view of the low temperature expansion, we are interested in 
the probability measure Pd-ising, on *p := UeeiJi(S;Z/2Z) ^(G) x V € (G), defined by, for every 

(-Fred) -Fblue) £ V : 

, , c (n ee p rcd e- 2 ^)(n ee p bluc e- 2 ^ 

-P > d-Ising(F > r ed) -Fblue) — 



-Ising 

(G*,J) 
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where C = (2]Q eg £. e J< =*) , and the partition function ^d-ising (G* , J) is given by: 

z d . lsmg (G*,j)= y: e c ( n e - 2Je *)( n ^) 

eeHi(E;Z/2Z) (P rcd ,P blue )eP e (G)xP E (G) eeP rcd eePbiue 

= E (^ing(^)) 2 - 

eeHi(S;Z/2Z) 

Given a pair (P rc d> Pblue) G ^Pj an d looking at the superimposition P rec j U Pbiuc on G, one 
defines two new edge configurations: 

• Mono(P re d, -Pblue): consisting of monochromatic edges of the superimposition P re d U 
Pblue) i- e - edges covered by exactly one of the polygon configuration; 

• Bi(P rc d, -Pbiue) : consisting of bichromatic edges of the superimposition, i.e. edges cov- 
ered by both polygon configurations. 

Edges which are not in the two configurations above are covered neither by Pblue nor by P re d- 
In sections 13.11 and 13.21 below, we characterize these two set of edges. 

3.1 Monochromatic edges 

Let e £ ifi(E; Z/2Z), and consider a pair of polygon configurations (P re d> Pblue) in V e (G) x 
V e {G). Then, it can be realized as four pairs of Ising spin configurations (±<r, ±<t')j each with 
defect type e, where coupling constants are negated along the representative e of e, chosen in 
Section 12.11 

Following Wilson [Wil li], to each of the four pairs of spin configurations, one assigns an 
XOR-spm configuration defined as follows: at every vertex, the XOR-spin is the product of 
the Ising-spins at that same vertex. 

Note that the four pairs of spin configurations yield two distinct XOR-spin configurations, 
one being obtained from the other by negating all spins. As a consequence, both XOR- 
spin configurations have the same polygon configuration separating clusters of ±1 spins, 
meaning that this polygon configuration is independent of the choice of (±<r, ±cr') realizing 
(-P re d) Pbhxe), let us denote it by XOR(P re d, -fWe)- Note also, that although the definition of 
a and a' depends on the particular choice of representative e, the XOR polygon configuration 
does not: it is defined intrinsically from (P re d) -Fblue)- 

Lemma 8. For every pair of polygon configurations (P re d, Pblue) ^ ^P; the monochromatic 
edge configuration Mono(P re d, Pblue) is exactly the XOR loop configuration XOR(P re< i, Pblue)- 
In particular, it is a polygon configuration ofV°(G). 

Proof. Fix a pair (-Pred; -Pblue) G of red and blue polygon configurations of G, that is 
(P red , Pblue) G V e {G) x V e (G) for some e G fli(£;Z/2Z). Let {a, a') be one of the four pairs 
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of spin configurations whose low temperature expansion is (P re d) Pblue)- We need to show 
that, for every edge e of G, e is monochromatic, if and only if XOR-spins at vertices u* , v* 
of the dual edge e* are distinct. 

Suppose first that e does not belong to e. Then by the correspondence between contours and 
spins given by the low temperature expansion, the edge e is monochromatic, if and only if 

o~ u * 7^ o~v* and a' u , = a v * , when e is only red, 
a u * = a v * and a' u * ^ a' v * , when e is only blue. 

If now e belongs to e, then for the same reasons, e is monochromatic if and only if one of the 
two above conditions holds, with colors exchanged. In all cases, the edge e is monochromatic 
if and only if XOR spins at vertices u* and v* are distinct. 

Being the boundary of some domain, the set of monochromatic edges must be a polygon 
configuration of G, with homology class 0. □ 



3.2 Bichromatic edge configurations 

Before describing features of bichromatic edge configurations, we recall some general facts. A 
polygon configuration P of the graph G separates the surface £ into np connected components 
Si,-- - , £ np , where np > 1. For every i G {1, • • ■ ,np}, £j is a surface of genus gi with 
boundary (9£j. The boundary is either empty or consists of cycles of £. 

As in Section [21 G^v denotes the subgraph of G, whose vertex set Vs ! is V n £j, and whose 
edge set E^. consists of edges of E joining vertices of Vz t , from which edges on the boundary 
dTii are removed. The dual graph is denoted by G^ . . 

Recall that Pi(£j, <9£j; Z/2Z) denotes the first homology group of Ej relative to its boundary. 
Consider the morphism 11^, from i^i(E;Z/2Z) to i?i(Ej, 9£j; Z/2Z) defined as follows: for 
every e G .ffi(£; Z/2Z), ITsj(e) is the homology class in Hi(Y<i, <9Ej; Z/2Z) of the restriction 
of any representative e of e to Ej, see Appendix IA.6I for details. 

The following lemma characterizes bichromatic edge configurations. 

Lemma 9. Fix e G f/i(E; Z/2Z), and let P G V°(G) be a polygon configuration, separating 
the surface E into connected components Ei, • • • , S np . 

• If there exists a pair of polygon configurations (P re d, Pblue) G T ,e (G) x V € (G) such that 
Mono(P re d, Pbiue) = Pj then, for every i G {1, • ■ ■ ,np}, the restriction of bichromatic 
edges to G^v is the low temperature expansion of an Ising configuration on G^,, with 
coupling constants (2J e *), and defect condition IIx;. (e). As a consequence, it is a polygon 
configuration in "P nE ^ e ^(GsJ- 

• Given, for every i G {1, • • • , np}, a polygon configuration Pi G V Us ^ e \Gj: i ) ; then, there 
are 2 np ~ l pairs (P re d, Pblue) £ V t {G) x V e (G) such that Mono(P re ^, Pbiue) = P and such 
that, for every i G {1, • • • ,np}, the restriction of bichromatic edges to G^ i is Pj. 
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Proof. 



• Suppose that there exists a pair of polygon configurations (P re d> Pbluc) of V e (G) x V e (G) 
such that Mono(i- > rc d, -Pblue) = P- Then, for every i G {l,--- , rip}, the restriction 
of bichromatic edges to G-£ { exactly consists in the restriction to of one of two 
original polygon configurations. Since this polygon configuration has homology e in E, 
the homology class in i?i(Ej, <9Ej; Z/2Z) of the restriction to Ej is Ils^e) by definition. 
As a consequence, the bichromatic edge configuration on Ej is a polygon configuration 
of P nE <^(G Ei ). Moreover, since all edges in the bichromatic configuration are present 
twice, and since the weight of pairs of polygon configurations is the product of the edge- 
weights contained in the pair of configurations, the effective weight of a bichromatic 
edge e is squared and becomes: 



which corresponds to a doubling of the coupling constants. 

• There are two spin configurations, denoted by ±£, whose low temperature expansion is 
P. Suppose that there exists a pair of spin configurations (a, a 1 ) whose low temperature 
expansion has P as monochromatic edges, then era' = ±£. Let us assume aa' = £, the 
argument being similar in the other case, this has the effect of adding a global factor 
2 when speaking of spin configurations. The relation a a' = £ implies that there is 
freedom of choice for exactly one spin configuration, say cr, the other being determined 
by their product £. 

Consider a connected component Ej, and a polygon configuration Pj G P nE ^ e ^(Gs 4 )- 
We want Pj to consist of doubled edges, so that in particular, it must contain all red 
edges. There are thus two choices for the first spin configuration of G s ., denoted by 
±cj 1 . This holds for every i G {1, ■ ■ ■ ,np} and thus defines 2 np spin configurations 
(iff 1 , • • • , ±a np ) of G*. Recall that in each of the 2 np cases, the second spin configu- 
ration is determined by the condition aa' = £. Since on each connected component Ej, 
£ is identically equal to ±1, we deduce that (a') 1 = ±cr\ As a consequence, the low 
temperature expansion of a' exactly consists of edges of Pj, i.e. Pi consists of red and 
blue edges. Summarizing, there are 2 • 2 np pairs of spin configurations, or 2 np_1 pairs 
of polygon configurations (P rc d, Pbiuc); such that monochromatic edges are those of P 
and bichromatic edges those of Pj, i G {1, • • • , rep}. Note that by construction (choice 
of a l, s), each polygon configuration P rc d, Pbluc is in V e (G). 



Let P 6 V°(G) be a monochromatic edge configuration, and let e G Pi(S, Z/2Z). Denote by 
WJ_ Isi [Mono = P] the contribution of the set 




□ 



{(Ped, Pbluc) G P £ (G) X P e (G) : Mono(P red , Pbluc) = P} 
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to the partition function (^i s i ng (<^)) 2 , and by 

W d . Ising [Mono = P} = ]T W d _ Ising [Mono = P) . 

eGHl(S;Z/2Z) 



By the low temperature expansion of the Ising partition function, the weight of each poly- 
gon configuration P-ed, Pblue is the product of edge- weights contained in the configuration. 
As a consequence, the contribution of (P re d> Pblue) can be decomposed as a product over 
monochromatic edges, and bichromatic edges of each of the components. Using Lemmas [8] 
and El this yields 

Proposition 10. For every P E V°(G) and every e E Pi(£, Z/2Z), 

np 

W d _ Ising [Mono = P] = 2- 1 C(H e- 2 ^)(j[2Z^ ie \G^2J)), (6) 

eSP 8=1 

where C = (2 Hee-E e<7e * ) 2 ■ Moreover, the double Ising partition function can be rewritten as: 

Zd-Ising(<7) = Wd-Ising[MonO = P], 

and the probability measure Pd-ising induces a probability measure on polygon configurations 
ofV°(G), given by: 

VP G P°(G), Fd- Ismg [Mono = P] = Wd - I - g[M °?r P] - (7) 



4 Mixed contour expansion 

In |Nie84j , Nienhuis rewrites the partition function of the Ashkin- Teller model on the square 
lattice as a statistical sum over polygon families on G and G* which do not intersect. We 
apply the same approach for the double Ising model on E but some care is required to keep 
track of the homology class of the polygon configurations involved. 

We fix e E fli(£;Z/2Z) and a polygon configuration P E V°(G). In Proposition [TTl we 

apply the low/high temperature duality to each of the terms Z L ^ e \G^ i , 2 J) involved in 
the expression of WJ_ Ising [Mono = P] of Equation ([6]). This has the effect of transforming 
bichromatic polygon configurations of G^. into polygon configurations of G^ ■ Then, in 
Proposition 1131 we sum over e E Pi(X; Z/2Z), and show that the outcome simplifies to a sum 
over polygon configurations of the dual graph, having homology class in Hi(J2; Z/2Z), and 
not intersecting P. 



14 



Proposition 11. For every e G Pu(£; Z/2Z) and every monochromatic polygon configuration 
P G V°(G), the contribution W,!j_ Ising [Mono = P] to (Zf sing ( J)) 2 , is equal to: 

cQiif^yni e (-i)<^ w ' e (n 

e£P i=l r<6Hi(S <) Z/2E) P/ST^G* .) e *GP* 

w/iere ; C = 2l v '*l+ 1 (n ee£ cosh(2J e . )) . 

Proof. The expression for W^_ Ising [Mono = P] of Equation flSJ) can be rewritten as: 

np 

WI_ Ising [Mono = P] = 2"-- 1 C(n e- 2J -) (II 2J)) . 

eeP 1=1 

For every i G {1, • • • ,np}, the contribution, Z^*^ {Gs i7 2J) is the low temperature expan- 
sion of an Ising model on vertices of VS. with coupling constants 2J e * and defect condition 
Ils^e). By using the relation between Kramers and Wannier's low and high temperature 
expansions of (|5|), it can be expressed as a product 

Z^ (e) (G Si ,2J) = A x £ (-l)^l n ^^Z 5 V(G Sl ,2J), 

r l ePi(£;;Z/2Z) 

where 

ilVS-l- 1 TT cosn ( 2 ^e*) 



A = 2 1 ^'- 1 n 

eeP E , 



e 2J e . 



Let us first compute the part which is independent of e. Observing that the collection of sets 
of dual vertices (V^.)™^ is a partition of V*, one writes: 



2^-^(11 e- 2Je *)(n A) =2^- 1 2 2 (ne 2 ^)([]e- 2JE *)2l v '*l-^( J] ^^). 

eeP i=l eeP eeP eeE ^i 

Noticing that the collection of edges in the Ej's is exactly the set of edges of G not in P, we 
have: 

up 

2 n r-lc(Y[e- 2J e*)(f[Ai) =2l v "l+ 1 ([]e 2J ^)(ne- 2Je *)( J] e" 2J <- cosh(2J e *)) 

eeP i=l eSE eSP e£P\P 

= 2l y *l +1 ( JJ cosh(2J e » )) (J] cosh(2 J e , )- x ) . 

eSE eeP 

Define the constant C = 2' v '*l +1 (]^[ eg£; cosh(2 J e *)) , then by definition of Z^ T (Gs i7 2J), one 
deduces that the contribution of Wj_ Ising [Mono = P] is equal to: 

np 

c'tn cosher 1 ) n[ E (-i) (rI|n ^ {e)) E (II tanh (2JeO) 

eeP i=l r'ePi(S I5 Z/2Z) P.*eP Ti (G|;.) e * eP i* 
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The proof of Proposition [IT] is concluded by observing that 

2e~ 2Je * 1 — e~ 4J - 

cosh(2J e *)~ = _ . , and tanh(2J f 



1 + e" 4J E * ' y e ' 1 + e~ 4J e* ' 

□ 

In order to have an explicit expression for the contribution of P to the partition function 
Zd-ising( J), we need to sum the quantities W|_ Ising [Mono = P] of Proposition [TT1 over e S 
i?i(S;Z/2Z). This is the object of the next proposition. 

Proposition 12. For every polygon configuration P £ V°(G), the contribution of P to 

w d , sing[ Mono = P ]= Cl (n rf^) e ( n • 

e£P {P*eP°(G*):PnP*=0} e*GP* 

w/iere Ci = 2l y *l+ 2 s +1 (]l ee£ cosh(2J e »)) . 

Proof. To simplify notations, let us write the product of weights of edges in polygon config- 
urations as follows: 

e£P e*eP* 

The contribution Wd-ising [Mono = P] can be rewritten as: 

n P 

W d . Isi n g [Mono = P]=C'e(P) E Il[ E (-l) (Tl|n ^ (£)) E *(^) • 

eePi(E;Z/2Z) «=1 r i 6ifi(E. i ,Z/2Z) P*^V T UG* ) 

Expanding the product over i 6 {1, • • ■ ,np} and exchanging the summation over e and 
(r 1 , • • • , r np ), one obtains that Wd-ising [Mono = P] is equal to: 

np 

CG(P) £ n[( e ®*(^))( e (-1)^-^1^^)". 

(^.•'■/" p )6n"i'iffiP.//2Z) i = 1 P*eV^(G^) eGPi(S;Z/2Z) 

Let us have a closer look at the last factor: YleeH 1 (T,-z/2Z)(~ l)^*=i( T '' ns ^ e ^. For every T = 
(r 1 , • • • , r np ) £ nr=i Z/2Z), consider the linear form: 

-0p: i?i(S;Z/2Z) — > Z/2Z 

_> Vr(e) = Er= P 1 (rin Ei ( e )), 

which once evaluated on the homology class e, gives the exponent of (—1) in the sum above. 
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Since (— 1)^t( £ ) equals 1 (respectively —1) if and only if e is (respectively is not) in ker(Y>7-), 
we can write: 

]T (-lfr(e) = #{e . ^ r(e) = 0} _ #{e . ^ r(e) = 1} 

eeHi(S;Z/2Z) 

= 2#ker(V> r )-##i(E;Z/2Z). (8) 

If -0T is a non constant form, its kernel is a hyperplane of -ffi(£; Z/2Z), thus having dimension 
2g—l and cardinality 2 2s_1 , implying that the sum over e of © is equal to 2.2 29_1 — 2 2g = 0. 
Therefore, the contribution of 7~ is non zero if and only if ipj- is identically zero, and in this 
case, the sum over e of (|8]) gives a factor: 

2#£Ti(E;Z/2Z) - Z/2Z) = #JETi(E;Z/2Z) = 2 2ff . 



Let us now determine which are the T's of (Jir=i -^i(Ei, Z/2Z)), such that ^7- is identically 
zero. 

First, notice that the evaluation of the intersection form (T*|Il£.(e)) on i?i(£j;Z/2Z) x 
Fi(S J ,aS i ;Z/2Z) is equal to (^(r^le) on JT a (S;Z/2Z) x #i(E;Z/2Z), where vrj is the pro- 
jection induced by the inclusion Ej C E, see Appendix IA.6I Indeed, take a representative r* 
of t* in Ej. Counting intersections with the restriction of a representative e of e to Ej is the 
same as counting intersections with the whole e, since t % is confined to Ej and thus has no 
intersection with e outside of Ej. Therefore, V't( c ) can be rewritten as: 

np np 
i=l i=l 

Since the intersection form is non degenerate on i?i(£; Z/2Z), ipj- is zero if and only if 
^"^-^(t 2 ) = G -ffi(£; Z/2Z), i.e. the homology class on the surface E of the whole 
polygon configuration P* = P£ U • • • U P* p is zero, i.e, P* G V°(G*). □ 

As an interesting corollary, we obtain a mixed contour expansion for the partition function 
of the double Ising model, and the corresponding expression for the double Ising probability 
measure of monochromatic polygon configurations which, by LemmaEl are the XOR polygon 
configurations, . 

Corollary 13. 

• The double Ising partition function can be rewritten as: 

z^ia-,j) = c E (IIiT^)(n 

{(P,P*)6P°(G)xP°(G*):PnP*=0} eeP e*eP* 

where Ci = 2l y *l +2 f+ 1 (J] eeB cosh(2J e »)) . 
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For every polygon configuration P £ V (G): 



-Ising[XOR — P] 



n 2e " 2J " ) y ( n 1 = e ~ 4Je * 

eGP {P*eP°(G*):PnP*=0} e*eP* 

z (n x + e -4j e , ) ( n i + e -4j e * ) 

{(P,P*)eP°(G)xpO( G *).p n p» = 0}. ee p e*6P* 



Note that one can see Kramers and Wannier's duality on this expression: the duality relation 
between coupling constant 

tanh J* = e~ 2J 

exchanges the expression for an edge of P and a dual edge of P* . 



5 Quadri-tilings and polygon configurations 

In the whole of this section, we let G be a graph embedded in a compact, orientable surface £ 
of genus g, and G* denote its dual graph. For the moment, we forget about the double Ising 
model. The goal of this section is to explicitly construct pairs of non-intersecting polygon 
configurations of G and G* , from a dimer model on a decorated, bipartite version G® of G, 
called quadri-tilings [dT07a]. 

This construction is done in two steps. The first step uses a mapping of Nienhuis |Nie84j . 
which constructs pairs of non-intersecting primal and dual polygon configurations, from 6- 
vertex configurations of the medial graph; this is the subject of Section [5TT1 The second step 
consists in using Wu-Lin/Dubedat's mapping |WL75t[Dubllb| from the 6-vertex model of the 
medial graph to the bipartite dimer model on the decorated graph G®. This is the subject 
of Section 15.21 

Using the above results and those of Section HI Theorem [18] proves that XOR loops of the 
double Ising model have the same law as primal polygon configurations of the bipartite dimer 
model. 



5.1 6-vertex model and polygon configurations 

The medial graph G M of the graph G is defined as follows. Vertices of G M correspond to 
edges of G. Two vertices of the medial graph are joined by an edge if the corresponding 
edges in the primal graph are incident. Observe that G M is also the medial graph of the dual 
graph G* , and that vertices of the medial graph all have degree four. Figure [3] represents the 
medial graph of a subset of 1? . 

A 6V- configuration or an ice-type configuration is an orientation of edges of G M , such that 
every vertex has exactly two incoming edges [Lie67j. An equivalent way of defining 6-vertex 
configurations uses edge configurations instead of orientations, as represented in Figure HI 
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Figure 3: The medial graph of Z 2 : plain lines represent Z 2 , dotted lines represent the dual 
graph Z 2 , and thick plain lines represent the medial graph (Z 2 ) M . Grey (resp. white) faces 
of the medial graph correspond to primal (resp. dual) vertices of the initial graph. 



This approach is more useful in our context, so that we define a 6-vertex configuration to 
be an edge configuration, such that around every vertex of G M , there is an even number of 
consecutive present edges. 




Figure 4: The six possible local configurations around a vertex v of the 6-vertex model, and 
their respective weights: arrow representation (top), and the even degree subgraph represen- 
tation (bottom). 

In order to make the 6-vertex model a model of statistical mechanics, weights are associated to 
local configurations around a vertex, and the probability of a 6-vertex configuration is taken 
to be proportional to the product of its local weights. In absence of external field, the weights 
of complementary local configurations are taken to be equal: there are thus three parameters 
for each vertex v of the medial graph G M , denoted by A v , B v and C v . Since, multiplying 
these three parameters by the same positive constant does not change the measure, we set 
C v = 1, see also Figure HI Let us denote by ZQ_ veTtex (G M , (A, B)) the partition function of 
this model. 

Mapping I [Nic84j. Consider the following combinatorial mapping from 6-vertex configu- 
rations to edge configurations of the primal and dual graph: whenever a vertex of G M has 
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two neighboring edges in the 6 vertex configuration, put the edge of G or G* separating the 
present and the absent edges, see Figure [5j The following lemma characterizes this mapping, 
see also [Nie84|. 



X 








M 
M 




Figure 5: Mapping on the local level (top), and on the global level (bottom). 



Lemma 14. 



• Mapping I associates to a 6-vertex configuration a pair of polygon configurations (P, P*), 
which do not intersect and such that the homology class of PUP* in i?i(S;Z/2Z) is 0. 

• Given a pair of polygon configurations (P,P*) as above, there are exactly two 6-vertex 
configurations which are mapped to (P, P* ) . 

Proof. A 6-vertex configuration consists of clusters of present/absent edges, see Figure \E\ 
The primal/dual edge configuration assigned by the mapping consists of the boundary of 
those clusters. It is thus a polygon configuration with homology class. A primal and the 
corresponding dual edge cannot intersect, since that would correspond to a configuration on 
G M where around a vertex there is alternatively one edge present, one edge absent, then again 
one present, one absent, which is a forbidden local configuration for the 6-vertex model. 

Conversely suppose we are given a pair of polygon configurations (P, P*) as above. The fact 
that the homology class of PU P* in Hi(Y,; Z/2Z) is 0, exactly means that it is the boundary 
of domains that can be painted alternatively in two colors consistently. Put all edges of G M 
in domains of one color, and remove all edges in domains of the other color. In this way, one 
exactly obtains two valid configurations of the 6-vertex model, one being the complement of 
the other, depending on which color is used to represent present edges, see also Figure [5j □ 

Let us denote by V°(G UG*) the set of pairs (P, P*) of polygon configurations of G and G* 
respectively, such that the union PU P* has homology class in .ffi(£; Z/2Z). 
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Recall that to every edge e of the primal graph (resp. e* of the dual graph), corresponds a 
vertex v of the medial graph, which we denote by v(e) (resp. v(e*)). Using this fact, one can 
naturally define a weight function on edges of G and G* : 

Ve G G, a e :=A v ( e y, Ve* G G* , b e * := B v ( e *y 

With this choice of weights and using Lemma [U we obtain the following. 

Lemma 15. 

Z 6 . vcrtcx (G M ,(A,B)) = 2 U a - II 6 e*- 

{(P,P*)&V°{GUG*): PnP*=0} eeP e*eP* 

5.2 Quadri-tilings and 6- vertex model 

Let us define yet another graph built from the graph G. The quadri-tiling graph of G, denoted 
by G®, is the decorated graph obtained from G M by replacing every vertex by a decoration 
which is a quadrangle. The graph G® is bipartite and can be drawn on the same surface 
as G. Edges shared by G® and G M are referred to as external edges, and those inside the 
decorations as internal. 

W' 

h* b e * 

,M, 

Figure 6: Decoration of a vertex of the medial graph G M (left) to obtain a piece of the 
quadri-tiling graph (right). 

A dimer configuration of G®, also known as a perfect matching is a spanning subgraph of 
G® where every vertex has degree exactly one. Let us denote by M(G®) the set of dimer 
configurations of the graph G®. In a particular decoration of G®, a dimer configuration of 
G looks like one of the seven possibilities represented in Figure [7] (top). 

Assigning positive weights (w e ) ef z E Q to edges of G®, the dimer Boltzmann measure, denoted 
Pq, is defined by: 

VMeM(G% P Q (M)= 7 Ue %™ e r , 

where Z^ mer {G ( ° i , w) = X^AfeM(GQ) IleGM w ^ ls * ne dimer partition function. This defines a 
model of statistical mechanics, called the dimer model on Cr*. 

Mapping II [WL75, Dubllb]. Requiring exterior edges to match yields a mapping from 
dimer configurations of G® to 6- vertex configurations of G M , see Figure [71 This mapping 
between local configurations is one to one except in the empty edge case where this mapping 
is two-to-one. 
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Figure 7: The local configurations of a quadri-tiling and the corresponding local 6-vertex 
configurations. 



We now choose weights of edges in a specific way. Let A and B be positive functions on 
vertices of the medial graph G M , defining weights of local 6-vertex configuration; and let a 
and b be the induced weight functions on edges of G and G* respectively. The weight function 



w 



a.b 



on edges of G® is defined as follows, see also Figure [6) 

1 if e is an external edge 
w a e ,h = < a s if e is an interior edge, parallel to a primal edge e 
be* if e is an interior edge, parallel to a dual edge e*. 

Let us denote by Zdi mcr (G^, (A, B)) the corresponding partition function. From now on, we 
suppose that local weights of the 6-vertex model satisfy the relation: 

Vu G 1/(G M ), A 2 v + B 2 v = l. 

This implies that, G V(G M ), A v = Av 2a"b Cv = ^> ^- e - the model is free fermionic. 
Using the weight functions a and b, the free fermionic condition can be rewritten as: 



Ve G E, a 2 e + b% 



1. 



(9) 



With this choice of weight, we thus obtain the following. 



Lemma 16. When local weights of the 6-vertex model satisfy the free fermionic relation, 
Mapping II is weight preserving and: 

Zdimer(G' Q , (A B)) = Z 6 _ vortex (G' M , (A, B)). 

Let us now choose weights of edges of G^ to depend on coupling constants (J e *) of the double 
Ising model: 



Ve G E, a e 



2e 
1 + e- 



-2J e » 
-4J„ 



and Ve* G E* 



, b e * 



1 + e- 



-4J e » 
^4jT 



(10) 
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Then, it is straightforward to check that a and b satisfy the free fermionic condition 

n o , 2 4e" 4J e* +(l-e- 4J E *) 2 
Vee£, a* + & = (1 + ^_ 4j£ , )2 =1. 



It should be noted that in the more general case of the Ashkin- Teller model, when spins of 
the two Ising models interact, the mapping with the 6- vertex model still holds |Mie84llSal87] . 
but the model is not a free fermion any more. 

Since weights a and b depend on coupling constants (J e *), we set J as argument for the 
corresponding dimer model partition function. 

Recall that the mixed contour expansion of the double Ising partition function, see Corol- 
lary [13l involves pairs of non intersecting primal and dual polygon configurations, each of 
which has homology class in Z/2Z). We want to take the same restriction here. 

Consider a dimer configuration M of the graph Gr*, then Mapping II assigns to M a 6- 
vertex configuration, and Mapping I assigns to this 6-vertex configuration a pair of non- 
intersecting polygon configuration of V°(G U G*). Let us denote this pair by Poly(M) = 
(Pol yi (M),Poly 2 (M)). 

We restrict ourselves to dimer configurations M such that Poly 1 (M) has homology class, 
and denote by Ai°(G^) this set. Note that since the superimposition Poly 1 (M) UPoly 2 (M) 
has homology class, this automatically implies that Poly 2 (M) also has homology. Let 
Pq be the corresponding dimer Boltzmann-measure, and ^dimer(^^' -0 be ^ ne corresponding 
partition function. 

Let (P,P*) G V°(G) x V°(G*) such that P n P* = 0. Denote by W Q [Poly = (P,P*)], the 
contribution of the set: 

{M G M(G Q ) : Poly(M) = (P,P*)} C M°(G), 

to the partition function Zj imer (G^, J). Then, as a consequence of Lemmas [T5l and [T6l we 
have: 

Proposition 17. When weights assigned to edges of the graph G® are chosen as in Equation 
<HDK, we have for all (P,P*) G P°(G) x P°{G*), such that P n P* = $: 

«^-)Mni&(n}^). 

eeP e*eP* 

Moreover, the dimer model partition function can be written as: 

j) - 2 E ttlif^XlIi^) 

{(P,P*)e-p°(G)x-p°(G*):PnP*=0} eeP e *eP* 
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and the probability measure Pq induces a probability measure on polygon configurations of 
V°(G), given by: 



»[Polyi = P] 



(ILeP x+ e -« e * ) „ (i Ic ••. /' |+, 

2e~ 2J e* \f-TT I < 



S (ILeP l+ e - 4 Je* ) ( rie*eP* l+e- 4J e 

{(P,P*)£V°(G)xP°(G*y. PnP*=V>} y ' x 



Combining Corollary 1131 and Proposition 1171 yields the following. 
Theorem 18. 

• The double Ising partition function and the dimer model partition function are equal up 
to an explicit constant: 

Z d .i smg (G*,J) = 2^l +2 s(n cosh(2J e ,))^ d ° im c r (G Q , J). 

eeP 

• XOR-polygon configurations of the double Ising model on G* have the same law as Polyx 
configurations of the corresponding dimer model on the bipartite graph G® : 

VP G V\G), P d _ Is mg[XOR = P] = F° Q [Pol yi = P]. 

Note that the first part of Theorem [TH] can also be deduced from the results of [Dub lib] and 
[CDT2] . 

Since the quadri-tiling model is a bipartite dimer model, it can be studied in great details 
the tools of Kasteleyn theory, which we recall some elements of in the next section. These 
tools can thus be used to study the distribution of the XOR Ising configurations. 



5.3 Kasteleyn theory 

We now recall some elements of the Kasteleyn theory for bipartite dimer models and apply 
it to the dimer model on £r*. This simplified version for bipartite graphs of the more gen- 
eral theory developed by Kasteleyn |Kas67| is due to Percus [Pe r69j . The main tool is the 
Kasteleyn matrix, defined as follows in the bipartite case: 

• rows (resp. columns) are indexed by white (resp. black) vertices; 

• the absolute value of an entry is if the corresponding white and black vertices are not 
adjacent, and is the dimer weight of the edge formed by these vertices when they are 
adjacent; 

• signs of the entries are chosen in such a way that around all (bounded) faces, the 
number of minus signs around a face has the same parity as half of the degree of the 
face, minus 1. 
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If the graph is planar, then the partition function of the model is, up to a global sign, the 
determinant of the Kasteleyn matrix. Indeed, when expanding the determinant of K as a sum 
over permutations, the only non zero terms are those corresponding to dimer configurations, 
and their absolute value is the correct weight. The third condition about signs is here to 
compensate the signatures of the permutations, so that all the terms exactly have the same 
sign. Kasteleyn showed [Kas61[ IKas67] that such a choice of signs exists (and is essentially 
unique): changing the signs of all edges around a particular vertex still yields a choice of 
signs satisfying the third condition, and one can pass from one valid choice to another by a 
succession of such operations. In order to get the right global sign, one can chose a reference 
dimer configuration Mq, giving a bijection between white and black vertices, agree that the 
order chosen for rows and columns of K is compatible with this bijection, and chose signs so 
that all entries of K corresponding to dimers of the reference configuration have sign +. 

When superimposing a dimer configuration M of with the reference dimer configuration 
Mo, each vertex of G^ is incident to exactly two edges of the superimposition, so that we 
obtain a family of non- intersecting loops and doubled edges covering all vertices of £r*. 

If the graph is embedded in a surface £ of genus g > 0, loops of the superimposition also 
live on £ and may have non-trivial homology in -ffi(£; Z/2Z). For the purpose of this paper, 
it is convenient to fix as reference dimer configuration Mq the one covering all internal edges 
of decorations parallel to dual edges of G* . Then we have the following useful lemma. 

Lemma 19. Let M be a dimer configuration ofG®. Then, the homology classes of MqVJ M 
and Poly 1 (M) m-ffi(£;Z/2Z) are equal. 

Proof. In order to prove that the homology classes are the same, it is sufficient to prove that 
they give the same result when computing the intersection form against any homology class 
r G Hi(£;Z/2%). 

Let us fix such a class r G Z/2Z). Let r be a representative of the class r, realized as a 

path on G* . We now show that the parity of the number of intersections between Poly 1 (M) 
and r, is equal to the number of intersection between MU Mq and r. All intersections occur 
in the interior of dual edge used by r. 

Fix e* a dual edge used by r. From the mappings above, the edge e belongs to Poly 1 (M) if 
and only if the number of interior edges parallel to e in the corresponding rhombus covered 
by dimers in M is odd (see Figures [5] and [7|) . Since edges of Mq are parallel to e* , the parity 
of the number of intersections with e* will be the same for Poly 1 (M) as for M U Mq. Since 
this holds for every dual edge belonging to r, it holds for r. Therefore, the homology classes 
for Poly 1 (M) and M U Mq are the same. 

□ 

Kasteleyn theory of graphs embedded in surface of genus g > is more complicated than for 
planar graphs. There is a topological obstruction for the existence of a sign distribution on 
edges giving every dimer configuration a + sign in the determinant expansion of a Kasteleyn 
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matrix. There still exist choices of signs satisfying the third condition for a Kasteleyn matrix, 
but there is not just one, as in the planar case, but 2 2g classes of choices of signs. Therefore 
we can define 2 2g non-equivalent Kasteleyn matrices. 

From one of them, denoted by , one can construct the others non-equivalent matrices de- 
noted by e G ffi(S; Z/2Z), by multiplying the sign of all edges crossing a representative 
of e by —1. 

In the expansion of the determinant of each of the 2 2g matrices, there are terms with different 
signs: The sign of a dimer configuration M in the expansion depends on e and the homology 
class in Z/2Z of the loops in the superimposition of M and of the reference dimer configuration 
Mq. With our choice of reference dimer configuration, it thus depends on the homology class 
of primal contours of Poly 1 (M). The determinant of each Kasteleyn matrix has thus the 
following form 

det V s Z {a) 

aefli(E;Z/ 2Z ) 

where is the partition function for dimer configurations such that Poly^^ has homology 

class a and s a>e is a sign depending only on e and a. 

It turns out that the linear relations between det K ^ and ^mer can ^ e inverted explicitly and 
each ^imer * s a nnear combination of the 2 29 determinants of Kasteleyn matrices [DZM + 96[ 
IGL991 ITesOOl ICR071 ICR08] . This is in particular true for ^- mer . 

Classically, when studying dimer models one is interested in the full partition function 

7 — ST 7^ 

^dimer — ^dimer' 

aeffi(£;Z/2Z) 

which is also a linear combination of the determinants of Kasteleyn matrices. 

Many properties of the dimer model with the full partition function are deduced from the 
study of Kasteleyn matrices. As a consequence, most of them can be readily obtained also 
for the restricted model where the homology class of configurations is fixed, since it is just a 
matter of considering another linear combination of determinants. 

This discussion of the relation of signs and homology considerations and the selection of 
particular homology classes via adapted linear combinations is not specific to the quadri- 
tiling model, but applies as is to any bipartite dimer. It applies in fact to any dimer model, 
where the more general Kasteleyn theory applies, using full Kasteleyn matrices with rows 
and columns indexed by all vertices of the graph, and Pfafnans instead of determinants. We 
refer the reader to [CR07| for an intrinsic geometric interpretation of coefficients in the linear 
combination in this general context. The low temperature polygon configurations of the Ising 
model can be mapped via Fisher's correspondence [Fis66] to a non-bipartite dimer model. 
Restricting the homology class of the polygon configurations can also be obtained on the 
dimer side with an appropriate linear combination of Pfafnans of Kasteleyn matrices. 
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6 The double Ising model at criticality on the whole plane 



After having discussed in much generality the case of finite surfaces of genus g, we now want 
to consider the case of infinite planar graphs. From now on, we restrict ourselves to a special 
kind of graphs, the so-called isoradial graphs, with specific values of the coupling constants 
for the Ising model. 



6.1 Isoradial graphs 

Definition 6.1. An isoradial graph [Ken02, KS05] is a planar graph with a proper embedding 
having the property that every face is inscribed in a circle of fixed radius, which can be taken 
equal to 1. 

The regular square, triangular and hexagonal lattices with their standard embedding are 
isoradial. A fancier example is given in Figure [8] (left). 




Figure 8: Left: an example of isoradial graph. Middle: the corresponding rhombus graph. 
Right: the half-rhombus angle 9 e associated to an edge e. 

The center of the circumscribing circle of a face can be identified with the corresponding dual 
vertex, implying that the dual of an isoradial graph is also isoradial. 

The dual of the medial graph G M , called the diamond graph of G and denoted by G , has as 
set of vertices the union of those of G and G* . There is an edge between v £ V and /* G V* 
if and only if v is on the boundary of the face / corresponding to the dual vertex /*, see 
Figured] (middle) for an example. Faces of G° are rhombi with edge length 1, diagonals of 
which correspond to an edge of G and its dual edge. To each edge e of G we can therefore 
associate a geometric angle 9 e G (0, 7r/2), which is the half-angle of the rhombus containing 
e, measured between e and the edge of the rhombus, see Figure [5] (right). The family (9 e ) 
encodes the geometry of the embedding of the isoradial graph. 

Note that in general, an isoradial graph has many isoradial embeddings, each of which gives 
a new collection of angles. 
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6.2 Statistical mechanics on isoradial graphs 

When defining a statistical mechanical model on an isoradial graph, it is natural to relate 
statistical weights to the geometry of the embedding, and thus to chose the parameters 
attached to an edge (coupling constants for Ising, probability to be open for percolation, 
weight of an edge for dimer models, conductances for spanning trees or random walk,. . . ) 
to be functions of the half-angle of that edge. For the Ising model, discrete integrability 
considerations (invariance under star-triangle transformations) and a self-duality argument 
yields the following expression for the interacting constants: 

J e = J{9 e ) = - log — . 

2 \ COS# e J 

This expression for J e for 9 e = j,^ and ^ coincides with the critical value of the square, 
hexagonal and triangular lattices respectively. The Ising model with these coupling constants 
and is known to exhibit a critical behavior [BDT10, BD T1H ICS09| . We therefore refer to 
these values as critical coupling constants. 

We suppose that the isoradial graph G is infinite, in the sense that the union of all rhombi 
of the diamond graph of G covers the whole plane. 

Consider the corresponding bipartite dimer model on the infinite decorated quadri-tiling 
graph G®. Then, the correspondence for the weights described in Section 15.2} Equation ([TP]) 
yields in this particular context: 

Ve 6 E, a e = a(0 e ) = cos6* e , Ve* G E* , b e * = b(6 e *) = cos6 e * = s'm9 e . 

Notice that with a particular embedding of the decoration, namely when external edges have 
length 0, and a and b edges form rectangles joining the mid-points of edges of each rhombus, 
the quadri-tiling graph G® is itself an isoradial graph with rhombi of edge length ^, and that 
weights (up to a global multiplicative factor of |) are those introduced by Keynon |Ken02j 
to define critical dimer models on isoradial graphs. 

It occurs that for the Ising and dimer models on infinite isoradial graphs with critical weights 
indicated above, it is possible to construct Gibbs probability measures |dT07a[ IBDT11] . 
extending the Boltzmann probability measures in the DLR sense: conditional on the config- 
uration of the model outside a given bounded region of the graph, the probability measure 
of a configuration inside the region is given by the Boltzmann probability measure defined 
by the weights above (and the proper boundary conditions). These measures have the won- 
derful property of locality: the probability of a local event only depends on the geometry 
of a neighborhood of the region where the event takes place. This means in particular that 
changing the isoradial graph outside of this region does not affect the probability. 

We can therefore consider the critical Ising model (resp. dimer, in particular quadri-tilings 
models) on a general infinite isoradial graph, as being that particular Gibbs probability 
measures on Ising configurations (resp. dimers configurations) of that infinite graph. 
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We now work with a fixed infinite isoradial graph G. We denote by IPj^ n „ the measure on 
configurations of the critical Ising model on G* . By taking two independent copies of the 
critical Ising model on G* , we get the Gibbs measure for the critical double Ising model 
^Ssing = Sng ® ^Sng' f rom which XOR contours can be constructed, as in the finite case. 
We denote by Pq the Gibbs measure on dimer configurations of the infinite graph Cr*. 

6.3 Loops of the critical XOR Ising model on isoradial graphs 

It turns out that the identity in law between polygon configurations of the critical XOR Ising 
model on G* and those of the corresponding bipartite dimer model on G® remains true in 
the context of infinite isoradial graphs at criticality: 

Theorem 20. Let G be an infinite isoradial graph. The measure induced on polygon config- 
urations of the critical XOR Ising model on G* , and the measure induced on primal contours 
of the corresponding critical bipartite dimer model on G® have the same law: for any finite 
subset of edges £ = {e%, . . . , e n }, 

PSsing[£ C XOR] = F$[£ C PolyJ. (11) 

Proof. We first prove (jlip when G is a bi-periodic isoradial graphs. In that case, the quotient 
of G by the translation lattice A can be viewed as an isoradial graph on the torus. Quotienting 
G by sublattices nA, with n G N + gives graphs G n on larger and larger tori, as n increases. 

If we apply Theorem 1181 to the specific case of the double Ising model on a toroidal, isoradial 
graph G* with critical coupling constants, we know that on G n , XOR polygon configurations 
have the same law as primal polygon configurations of the corresponding bipartite dimer 
model on Gn- 

On the one hand, the graph G® being bipartite, we know [KOS06J that the sequence of 
Boltzmann measures Pq" of the dimer model on Gn , without any restriction on the homology 
of (primal) cycles, converges to the Gibbs measure Pq for quadri-tilings on G®, studied in 
[dT07a]. The same arguments can be used to show that the measure on even configurations 
of Gn also converges to the same limit Pq. 

On the other hand, the low temperature contours on G n for the Ising model on G* n are mapped 
to a (non-bipartite) dimer model via Fisher's correspondence [Fis66j. Using Kasteleyn's 
theory, the authors have proved in [BDT10], Theorem 6, that the law of the low temperature 
contours on G n converges to the Gibbs measure P^? on contours of G. In the proof of this 
theorem, the contours considered have no constraint on their homology class, but a careful 
analysis is carried out for each of the four terms appearing when expressing the measure in the 
Fisher dimer setting. Since imposing the homology class of contours, just as in Kasteleyn's 
theory for bipartite graphs, is just a matter of taking another linear combination of these 
four terms, the analysis goes through and the limiting measure exists and is the same as the 
non restricted case whatever the topological defect for the Ising model is. As a consequence, 
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the probability measure on the double critical Ising model converges weakly to the product 
measure P^ Ising . 

Since containing £ for XOR (resp. primal) polygon configurations is a local event, i.e. a finite 
union of cylinder events, for the double Ising model (resp. for the bipartite dimer model on 
G®), the equality between probabilities of £ , which holds on G n for every n, also holds in the 
limit: this implies that the law of the polygonal configuration of the critical XOR Ising model 
on the infinite isoradial graph G is the same as the law of the primal polygon configuration 
in the critical dimer model on Gr* 3 which ends the proof of the theorem in the case where G 
is bi-periodic. 

Furthermore both measures, for the critical Ising model on G* [BDTllJ and the critical dimer 
model on G® [Ken02[ TdT07a] . have the locality property: in both models, the probability of 
a local event only depend on the geometry of the embedding of G in a domain containing the 
spins or edges involved in this event. This is the crucial property required for the construction 
of the Gibbs measure of these models on non periodic infinite isoradial graphs [dT07aj. It 
implies that one can identify the probability of local events on a non periodic graph G with 
the one on periodic graphs with large fundamental domains coinciding with G on a ball 
sufficiently large to contain a neighborhood of the region of the graph involved in the event. 

As a consequence, probabilities still coincide in the non periodic case: the equality in law 
holds also when G is not necessarily bi-periodic. □ 

7 Height function on quadri-tilings 

Dimer configurations of the quadri-tiling graph G®, like all bipartite planar dimer models, 
can be interpreted as random surfaces, via a height function. It is the main ingredient to 
relate the previous results connecting XOR loops and dimers with Wilson's conjecture. 

7.1 Definition and properties of the height function 

Let us now recall the definition of height function, used in }dT07a) . A dimer configuration 
M of a planar bipartite graph can be interpreted as a unit flow o>m, flowing by 1 along each 
matched edge of M , from the white vertex to the black one. It is a function on edges having 
divergence +1 at each white vertex and —1 at each black vertex. Subtracting from ctjf 
another flow with the same divergence at every vertex, yields a divergence-free flow, whose 
dual is the differential of a function on faces of this graph. 

There is a natural candidate for this unit reference flow: since in a dimer configuration there 
is exactly one dimer incident to every vertex, the sum over all edges incident to any given 
vertex of the probability that this edge is covered by a dimer, is equal to 1. This means that 
the flow a>o, flowing by Pq (ej3 from the white vertex to the black one along each edge e of 
the graph, is a flow with divergence +1 (resp. —1) at every white (resp. black) vertex. 

1 The graph is isoradial and infinite, and the weights for the quadri-tilings are critical. So in this 
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The height function h on quadri-tilings is defined as follows. For every dimer configuration 
M of G^, h M is a function on faces of Cr*, such that for every pair of neighboring faces / 
and /' of sharing an edge e, with the additional property that when traversing e from / 
to /', the black vertex of e is on the left: 

h M (f')-h M (f) = a M (e)-a (e). 

When faces / and /' are not incident, choose a path / = /o, /i, • • ■ , f n = f in the dual graph 
joining / and /', then: 

n-1 

h M (f) - h M (f) = 5>"(/*+i) - h M (f i )). 
i=0 

This definition is consistent, i.e independent of the choice of path from / to /', because the 
flow Qj\j — ctQ is divergence free; it determines h up to a global additive constant, which 
can be fixed by saying that the height at a particular given face of G Q is 0. Faces of G Q arc 
split into three distinct subsets: vertices in the center of faces of G°, vertices of the primal 
graph G and vertices of the dual graph G* . We suppose for the sake of definiteness that the 
face where the height is fixed at corresponds to some particular vertex of G. 

Note that another choice of reference unit flow could have been the one coming from the 
reference dimer configuration Mq, where a white-to-black unit is flowing along all interior 
edges parallel to edges of G* . This would have produced a random height function whose 
restriction to the vertices of G (and to those of G*) would have coincided with that of this 
one (up to an additive constant). 

Denote by hy (resp. hy*) the restriction of h M to vertices of G (resp. to vertices of G*). 

The next lemma describes possible height changes between pairs of vertices of the primal 
(resp. dual) graph, incident in the primal (resp. dual) graph. To simplify the picture, we 
consider primal and dual vertices to be around a rhombus of the diamond graph, see Figure [9) 

Lemma 21. The following assertions are true: 

• the function hy (resp. h^* takes values in Z (resp. Z + \); 

• The increment of hy (resp. hy*) between two neighboring vertices of G (resp. of G* ) 
is —1, 0, or 1; 

• The increment ofhy (resp. hy*) is non zero if and only if the two vertices are separated 
by an edge o/Poly 2 (M) (resp. Poly 1 (M)j. 

particular context, we know |Ken02j that the probability of an edge is given by 9/n, where 9 is the half-angle 
of the rhombus containing that edge. 
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Figure 9: Height changes for the dimer model in a rhombus of the diamond graph. 



Proof. Let ei, e2 be two interior edges of C**, parallel to an edge of G, as in Figure [9) Then, 
the reference flow «o has the same value but opposite direction on these two edges. As a 
consequence, using the definition of the height function, 

A similar expression holds for /i M (v 2 )-h M (vi). This proves that the increment of h M between 
two neighboring vertices of G (resp. G*) is equal to —1, 0, or 1. Because of our convention 
for the base point, this implies that h M takes integer values on G. To see that h M takes 
half-integer values on G*, one just has to notice that the reference flow ao separating two 
vertices v (on G) and v* (on G*) which are neighbors on G° is = \ since the corresponding 
rhombus in the isoradial graph G® is flat. □ 

The level lines of hy (resp. hy*) are the set of closed contours on G* (on resp G) separating 
clusters of vertices of G (resp. of G*) where hy (resp. hy*) takes the same value. 

Returning to the definition of the pair of polygon configurations Poly(M) assigned to a 
quadri-tiling M, we immediately obtain the following: 

Lemma 22. Let M be a dimer configuration ofG®, then level lines of hy , respectively hy*, 
exactly correspond to the polygon configuration Poly 1 (M) ; respectively Poly 2 (M). 

Note that due to the fact that Poly 1 (M) and Poly 2 (M) do not cross, the increments of h M 
along two diagonals of a rhombus cannot be both non zero. As a consequence, on contour 
lines of hy , hy* is constant. 

Combining Lemma [22] with Theorem 1201 stating that monochromatic polygon configurations 
of the XOR Ising model have the same distribution as primal polygon configurations of dimer 
configurations of Cr*, we obtain one of the main theorem of this paper: 
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Theorem 23. Monochromatic polygon configurations of the critical XOR-Ising model have 
the same distribution as level lines of the restriction to primal vertices of the height function 
of dimer configurations of G® . 

7.2 Wilson's conjecture 

In [Willi j . Wilson presented extensive numerical simulations on loops of the critical XOR 
Ising model on the honeycomb lattice, on the base of which conjectured the following: 

Conjecture 1 (Wilson [Willi] ). The scaling limit of the family of loops of the critical XOR 
Ising model are the level lines of the Gaussian Free Field corresponding to levels that are odd 
multiples of 

The level lines of the Gaussian Free Field corresponding to levels that are odd multiples 
of A = form a CLE4 [SS09J. The contour lines of the XOR Ising models are thus 
conjectured to have the same limiting behavior as the CLE4, except that there are \/2 fewer 
loops in the XOR Ising picture. This conjecture is in agreement with predictions of conformal 
field theory jiRTTl IPSTT] . 

Theorem [23] can be seen as a proof of a version of Wilson's conjecture in a discrete setting, 
before passing to the scaling limit, and brings some elements for the complete proof of this 
conjecture. In particular, it explains the link with the Gaussian Free Field and the factor 
y/2, as we will explain now. 

For e > 0, denote by G^ the embedding of G® in the plane where rhombi of have side 
length e. For every dual vertex v in G^ , define v £ the vertex in Ge corresponding to the 
dual vertex v. 

The random height function h can be interpreted on G® as a random distribution, i.e. 
a continuous linear form on the set C^° C (R 2 ) of compactly supported smooth, zero mean 
functions, denoted by H £ : for every ip £ Cq° c (IR 2 ), 

H £ (<p) = ^ ™ea(v £ )h(v)^(v £ ), 

where area(t> £ ) = e 2 area(t>) is the area of the face of G® associated to v £ . 

In [dT07b] , the second author proved the following convergence result for the height function 
of the dimer model on G^: 

Theorem 24 ([dT07b]). As £ goes to 0, the height function on the critical quadri-tilings, as 
a random distribution, converges in law to times the Gaussian Free field. 

The result also holds for the restriction of h to G (resp. to G*) as soon as area(v) is replaced 
by the area of the corresponding face of G (resp. G*). 
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As contour lines of the restriction of h to G separate integer values, they can be understood as 
discrete level lines corresponding to half-integer values. Therefore, it is natural to expect that 
these contour lines converge to the contour lines of the limiting object, i.e to level lines for 
the Gaussian Free Field with levels (k + ^)y/n, fc£Z, which would prove Wilson's conjecture. 
Unfortunately, the result for the convergence result of the height function to the Gaussian 
Free Field is too weak to ensure convergence of contour lines. 

The convergence result in the paper [dT07b| applies not only to critical quadri-tilings, but to 
all bipartite planar dimer models on isoradial graphs with critical weights. It is conjectured 
that the family of loops obtained by superimposing two independent critical dimer configura- 
tions converges to CLE4. This is supported by the fact that each of the dimer configurations 
can be described by a height function, converging in the scaling limit to 1/y/ir times the 
Gaussian Free Field, the two fields being independent. Dimer loops are the half integer level 
lines of the difference, which by independence converges (in a weak sense) to y/2/ir times the 
Gaussian Free Field, and it is known that level lines (k + 1/2) of the Gaussian Free Field 
are a CLE4. 

Therefore, the factor y/2 in Wilson's conjecture corresponds to the fact that contours in the 
XOR Ising model have to do with contour lines of only one dimer height function, as opposed 
to two for dimer loops. 

As a consequence, 

Proposition 25. Conditional on a proof of the convergence of dimer loops to the CLE4, 
Wilson's conjecture is true. 

A Some elements of homology theory on surfaces 

Here are some general facts about homology theory on surfaces which are useful in the context 
of this paper. More details can be found in the reference [Ful951 [Mau96} IMas91] . We consider 
E to be a compact orientable surface of genus g with boundary <9E consisting of p components. 
The boundary maybe empty, in which case p = 0. 

We are interested in the first homology group Hi, and in the case where the target abelian 
group is Z/2Z. The other non trivial homology groups -Ho an d H2 are isomorphic to Z/2Z, 
when E is connected. 

A.l 1-chains and first homology group 

A 1-chain is a 1-dimensional submanifold of E. A 1-chain is a cycle if its boundary is empty. 
Among the cycles, there are those which are the boundary of a 2-dimensional submanifold 
of E. Those are simply called boundaries. The first homology group i?i(E;Z/2Z) is the free 
abelian group generated by cycles, modulo boundaries: two cycles 7 and 7' represent the 
same element in i?i(E;Z/2Z) if their union is the boundary of a 2-dimensional submanifold 
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of E. Note that since the target group is Z/2Z, we do not need to care about orientations of 
1-chains. 

A. 2 Relative homology 

If A is a closed subset of E, then one can also consider homology relative to A. As above, one 
defines a notion of cycle and boundary, relative to A this time: a relative cycle is a 1-chain 
whose boundary is in A. A relative boundary is a 1-chain in S for which there exists a 1-chain 
in A, such that the union of the two is a boundary in E. Then, the first homology group 
relative to A, denoted by i?i(E, A; Z/2Z), is the free abelian group generated by relative 
cycles, modulo relative boundaries. 

A particular example of interest is when A = 9E. Representatives of equivalence classes 
of i?i(E, <9E; Z/2Z) are finite unions of cycles and paths attached to components of the 
boundary. 

Note that when <9E is empty, the homology of E relative to its boundary coincides with the 
usual homology. 

A. 3 Explicit bases of homology 

The two homology groups i?i(E;Z/2Z) and i?i(S, <9E; Z/2Z) turn out to have the same 
dimension N = 2g + p — 1. For each of the groups, a basis can be explicitly given. Label 
the handles of E from 1 to g, and the p components of the boundary from Cq to C p _i. For 
iii(E;Z/2Z), choose N cycles (Aj)^ on E, as follows: 

• for i G {1,... ,<?}, take A 2 j-i and A 2 j to be winding around the i-th handle in two 
transverse directions, 

• for i E {1, . . . ,p — 1}, take A 29+ j to be wind around Cj, without crossing Ai, . . . , A 2g . 

Denote by Aj the homology class of Aj. Then, the collection (Aj)^ 1 is a basis of -Hi(E; Z/2Z): 
the homology class of any 1-chain on E is homologous to a sum of Aj's. The first homology 
group i?i(E;Z/2Z) is isomorphic to (Z/2Z) Ar : for every i 6 {1, • • • ,N}, the basis element Aj 
is mapped to the basis element of (Z/2Z) Ar consisting of O's and a 1 at position i. 

For #i(E,dE;Z/2Z), choose N cycles (7;)^ l on E, as follows: 

• for i G {1, ... ,5}, take 7 2i -i = A 2i , and -y 2 i = A 2i -i, 

• for i € {1, . . . ,p — 1}, take ~f2g+i to be a path from Co to Cj. 

Denote by 7^ the relative homology class of 7^. Then, the collection ( , ji)fL 1 is a basis for 
iJi(E, <9E; Z/2Z). The group ffi(E, <9E; Z/2Z) is also isomorphic to (Z/2Z) Ar . The two bases 
(Xi)fL 1 and (7i)^Li are dual to each other as explained in Appendix IA.5I 
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A. 4 Representatives of homology classes on graphs 



Consider a discretization of the surface E by a graph Gs, as in Section [2j The embedding 
of Gy, on T defines a notion of dual graph for Gs, denoted by G^. Then representatives of 
any homology class of i?i(E;Z/2Z) (resp. any relative homology class of H\(T, dT; Z/2Z)) 
can be realized as combinatorial paths on G* (resp. G). Figure [TOl provides an example of 
representatives of the bases (Xi)fL 1 and (74)^ defined in Section IA.3I 




A. 5 Intersection form 

There is a natural pairing between H\(T; Z/2Z) and H\(T, dT; Z/2Z), called the intersection 
form: 

: ffi(S;Z/2Z) x fl^E, <9E; Z/2Z) — ► Z/2Z, 

defined as follows. Let r £ JJj(E;Z/2Z) and e G H\(T, dT,; Z/2Z) be two homology classes. 
Take (generic) representatives r and e for two classes r and e respectively. Then (r|e) is 
defined as the parity of the number of intersections of r and e. This definition does not 
depend on the choice of representatives. 

In the explicit bases of Hi(T; Z/2Z) and Hi(T, dT; Z/2Z) chosen in Appendix lA.3[ the matrix 
of the intersection form is the identity. The pairing is thus non degenerate and defines an 
isomorphism between i?i(E; Z/2Z) and Hi(T, dT; Z/2Z). This is an explicit realization of 
the Poincare-Lefschetz duality []. 
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A. 6 Inclusion, excision and morphisms for homology 

Suppose that there exists a larger surface S containing S. Then a 1-chain in £ is in particular 
a chain in S, and a boundary in S is in particular a boundary in S. This means that the 
inclusion S C S induces a morphism 

7T S : iJi(S;Z/2Z) — »• Hi(£;Z/2Z). 

The inclusion also induces morphisms for relative homology groups: if the subset i C S 
is included in a subset B C S, then any relative chain (resp. cycle) in S relative to ^4 is 
in particular a relative chain (resp. cycle) in £ relative to B (just by forgetting what is in 
B\A). Therefore, this induces a morphism 

fli(E,A;Z/2Z) — ► ffi(E,5;Z/2Z), 

giving the homology class in £ relative to 5 of the restriction to £ \ S of any representative 
of an element of -Hi(£, A; Z/2Z). 

Moreover, the excision theorem states that if we cut out an open set U from both £ and 
A, the relative homology groups i?i(£, ^4; Z/2Z) and iii(£ \ £/, A \ J7; Z/2Z) are isomorphic. 
In particular, when J7 = £ c = £ \ S and ^4 = [7, then the excision theorem states that 
Fi(S,S c ;Z/2Z) and 9S; Z/2Z) are isomorphic. 

When composed with the previous morphism given by the inclusion E c C S, C S, it gives 
the following morphism 

n E : iTi(S;Z/2Z) — »• ffi(S,SE;Z/2Z). 
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